Critical currents and self organization in quantum Hall bilayers 
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We present a theory of the critical interlayer tunneling current in a disordered quantum Hall 
bilayer at total filling factor one, allowing for the effect of static vortices. In agreement with recent 
experiments [Phys. Rev. B 80, 165120 (2009); ibid. 78, 075302 (2008)], we find that this critical 
current is proportional to the sample area and is comparable in magnitude to observed values. This 
reflects the formation of a Bean critical state as a result of current injection at the boundary. We 
predict a crossover to a critical current proportional to the square-root of the area in smaller samples. 
We also predict a peak in the critical current as the electron density varies at fixed layer separation. 



In a quantum Hall bilayer at total Landau level filling 
Vt = lj Coulomb interactions induce a state with inter- 
layer phase coherence [l], Q , which can be understood as 
a Bose-Einstein condensate of interlayer excitons 0, 0] ■ 
The motion of excitons corresponds to counterflowing 
electrical currents in the layers, so that excitonic super- 
currents can give dissipationless electrical transport. In 
particular, a finite interlayer current I at negligible in- 
terlayer voltage V has been predicted @, @ and recently 
confirmed by four-terminal measurements 0, B- This 
behavior, which is a form of Josephson effect [9(, also 
appears as a zero-bias peak in differential conductivity 
measurements if 
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The low-voltage regime termi- 
nates when the interlayer current exceeds a critical value 
I c . Dissipation increases dramatically above I c . 

A key question has emerged over the nature of the 
critical current I c . Recent experiments have shown that 
it is proportional to the sample area [H, [l2|. Within 
the simplest model of a clean homogeneous bilayer, this 
scaling can only be explained by assuming that the tun- 
nel splitting is several orders of magnitude smaller than 
expected [13|. For realistic splittings, tunneling is es- 
timated to occur within a few microns of the contact 
[Eq. dU] so that I c should not depend on the sample 
length in the direction of the current 13- IH]. Edge tun- 
neling [l(| has also been proposed, predicting linear scal- 
ing with the sample length but not its area. 

In this paper, we present a theory which produces a 
critical current [Eq. (0] that is proportional to the area 
of the sample and of the correct order of magnitude, given 
reasonable estimates for the parameters. A disorder- 
induced lengthscalc, L^, emerges in our theory [Eq. (151)1 : 
this scale has no counterpart in the clean system jl3| . 
Our results are also consistent with the observed depen- 
dence of I c on the magnetic length and on an in-plane 
magnetic field. A key test of our theory is the prediction 
that I c should scale with the square-root of the sample 
size for samples smaller than 

The essential feature of our work is that we allow for 
static vortices in the exciton superfluid, which will be nu- 
cleated by strong charge disorder jl7H20i |. Such vortices 



play a crucial role for the critical current. They pin any 
injected supercurrents and sustain dissipationless states, 
in much the same way that disorder pins magnetic flux 
in superconductors [2l|-|23|. or charge in charge-density 
waves 24 1 . However, there is a significant difference in 



the bilayer: the depinning force comes from the injected 
charge current. Since these currents cannot penetrate 
the bulk of the quantum Hall state, the depinning force 
is applied only at the sample boundary. Given this ge- 
ometry, it is a surprising feature of our results that the 
critical current I c can scale with the sample area. 

We will first present numerical results showing that, in 
one dimension, currents injected at the boundary decay 
linearly in space in the disordered state (Fig. HJ) ■ This is 
analogous to the Bean critical state in a superconductor. 
We argue that these special critical states are generated 
by current injection from the boundary, and hence are 
selected in the bilayer. Based on this numerical work, 
we will then present a heuristic argument for the criti- 
cal current, which we generalize to two dimensions and 
compare with the experiments. 

Our starting point is the energy functional 
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for the low-energy modes of a bilayer containingpinned 
vortices. This form follows from the clean model |5| when 
the phase field of the pinned vortices 9° is subtracted out 
of the superfluid phase 6: <f> = 9 — 6 . The first term 
in Eq. (fTJ) is the superfluid stiffness while the second 
describes the interlayer tunneling. We assume that the 
vortex field 9° is disordered, with a correlation length £ 
which we estimate to be of the order of the distance to 
the dopants, dd ~ 200 nm in current samples [17j . The 
countcrflow supercurrent density above the ground state, 
jcF, and the interlayer tunneling current density, jt, are 
related to the phase field by: 
ep 
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A time-varying superfluid phase <f>(t) gives rise to an 
interlayer voltage difference V via the Josephson relation 
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V = hcb/e. Therefore, a state with a finite interlayer cur- 
rent at zero interlayer voltage is time- independent, corre- 
sponding to a local minimum of the energy ([1]) . To inves- 
tigate this possibility, we consider the dissipative model 
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-/9 s V 2 + tsin(0 + 6l ), 



(3) 



whose stationary solutions = V = are the local min- 
ima of Eq. ([TJ. The stationary equation is the conti- 
nuity equation stating that the loss of countcrflow cur- 
rent (first term) is accounted for by interlayer tunneling 
(second term). The dissipative dynamics of Eq. (|3|) is 
physically a resistive shunt due to interlayer quasiparti- 
cle tunneling. Additional dynamical terms are necessary 
in the finite voltage regime [T3, O US HHJ. However, it 
is the dissipation that determines the long-time limit at 
low voltages, so that here we may use Eq. (J3J). 

The boundary conditions for Eq. ([!]) come from the 
current flows through the sample [13l |. For definiteness, 
we consider a tunneling geometry in which a current It 
is injected into the top layer at one corner and removed 
from the bottom layer at the opposite corner. These 
current flows may be written as superpositions of layer- 
symmetric and layer-antisymmetric currents, 

/i„(out) = ! [(M) ±(1,-1)] 

where the components refer to the two layers. Thus, 
the tunneling experiment corresponds to a flow of layer- 
symmetric current, with equal countcrflow currents / = 
It/2 injected by both the electron source and drain. In 
the low-voltage regime, these countcrflow currents will 
be the supercurrents j'cf in Eq. ([2]). Since the symmetric 
component cannot penetrate the bulk, the supercurrents 
are injected at the boundary. As we shall see, the profile 
of their injection along the boundary is unimportant. 

The mechanisms leading to an extensive critical cur- 
rent for the bilayer can be seen most clearly from a one- 
dimensional model. In Fig. [TJ we show current and volt- 
age profiles obtained for the one-dimensional version of 
Eq. ([3]), on a lattice of L = 200 sites. These results 
are averaged over realizations of the disorder 6®, which is 
taken independently on different lattice sites. This corre- 
sponds to using the correlation length £ of the continuum 
model as the lattice spacing. The natural unit of current 
is then Iq = ep s /h£. For this illustration, we take the 
tunneling strength t£ 2 /p s = 0.4. In each realization, we 
start from an initial state in which the fa are random and 
independent, and equilibrate by integrating forwards in 
time with the boundary conditions d x fai = d x fa\L = 0. 
To model the current injection in a tunneling experiment, 
we then slowly increase the boundary conditions to the 
final values £d x <t>\i = —£,d x 4>\L = 1/1$. For the low- 
est three values of / used, the dynamics reach a time- 
independent solution, corresponding to the Joscphson 




200 



FIG. 1: Profile of counterflow current / (lines) and interlayer 
voltage (V oc <fi, vertical bars, non-zero in top curve only) from 
Q on a one-dimensional lattice with current injection at both 
ends, at a time w 10 4 /A after the current is switched on. The 
injected counterflow currents at each boundary are 3, 6, 9, 12/q 
for the four curves (la = ep s /h£/). The interlayer voltages 
vanish below a critical current I c . Here, 9 < I c /Io < 12. 
t£ 2 1 ' p s = 0.4, results averaged over 50 realizations. 



regime with vanishing interlayer voltages. For too large 
/, these time-independent solutions break down and the 
phase winds continuously at late times. This corresponds 
to the breakdown of the d.c. Josephson regime and the 
appearance of a state with finite interlayer voltages. 

The static states in F ig. [1] d iffer qualitatively from 
those of the clean model [lj, |26[, 8° = 0. In that case, 
there is a penetration depth for the injected current of 



Aj ~ V 'ps/t ■ 



(4) 



Since the phase angle is periodic this implies a maximum 
injected current density of d x <j> ~ 7r/Aj. In the disor- 
dered case shown, however, the injected current decays 
linearly close to the contacts, with a slope which is inde- 
pendent of the injected current. Thus, an increase in the 
injected current is accommodated by an increased cur- 
rent penetration into the sample. As can be seen from 
the figure, this process continues until the currents fill 
the entire sample. Beyond this point, further increases 
in current cannot be accommodated by coherent tunnel- 
ing and an interlayer voltage develops. Since the current 
decays linearly with a constant slope, the resulting I c in 
the one-dimensional model scales with the sample length. 

In the clean model, the breakdown of the stationary 
solutions can be understood in terms of the injection of 
phase solitons at the boundary [H, 27|, which propagate 
through the sample. Thus, the phase at any point varies 
in time, and the system develops an interlayer voltage 
by the a.c. Josephson effect. However, for the disordered 
system, the solitons may be pinned by disorder. We now 
develop a heuristic theory of such pinning, which agrees 
with our numerical work. 

We begin by recalling 24, 28| the form of the ground 



states of the random field XY model, Eq. ([T]) , in the weak 
disorder regime £ <C Aj relevant for the bilayer. In this 
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regime, the ground state consists of ferromagnetic do- 
mains with polarized phase. The key idea is that it is 
energetically costly to have phase twists at scales shorter 
than the size, Ld, of these domains. The energy cost for 
a phase twist that varies over the scale Ld is p s L®~ 2 in 
D dimensions. The typical tunneling energy of a polar- 
ized domain is obtained by summing random energies in 
the range ±t£ D for its (L d /^) D correlation areas, giving 
ti D {L d /i) D / 2 . Thus, for d < 4, the phase stiffness wins 
at short scales, and the domains reach a finite size where 
the two energies balance. This gives the Imry-Ma scale: 



La 
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(5) 



In this ground state of polarized domains, the average 
coarse-grained phase over a domain is chosen such that 
the tunneling energy H t of each domain is minimized. 
Since SH t /S(j)(r) is the tunneling current at position r, 
the total tunneling current over the domain vanishes. 

We now consider how the system changes in the pres- 
ence of an injected current. The injected counterfiow will 
cause the phase to twist away from its equilibrium value, 
leading to finite tunneling currents. We assume that the 
configuration remains smooth on the scale Ld, and hence 
average Eq. (|3]) over each domain. The second term in 
Eq. © becomes {t/L%)(L d /£) D / 2 f{4>) where 4> is the de- 
viation of the coarse-grained phase from its equilibrium 
value, and the range of /(</>) is typically [—1,1]. For a 
dissipationless state (0 = 0), the coarse-grained phase <j> 
should therefore obey: 



L 2 V 2 + f{$) = 



(6) 



The source term in Eq. © describes the loss of in- 
jected current due to tunneling in a domain. As discussed 
above, current injection induces counterfiow currents and 
hence phase twists. Since it is energetically costly to 
introduce phase twists in a domain, the domain at the 
boundary will respond by rotating uniformly, increasing 
its tunneling current, thereby reducing the counterfiow 
current. This process continues until the tunneling in 
the domain saturates, so that |/| ~ 1. The residual 
counterfiow currents will be transmitted further into the 
sample, causing the domains there to rotate in a similar 
way. Thus, we argue that forcing at a boundary leads to 
a self-organized critical state, in which the driven part of 
the system sits at the threshold |/| ~ 1. 

In one dimension, this argument means that Eq. ((6]) 
would give an average counterfiow current /3 s V</> that de- 
creases linearly from the boundary. To be more precise, 
it predicts a linear I/Iq in the saturated regions, with 
a slope — £ _1 (i£ 2 // 5 s) 4 ■ This is qualitatively consistent 
with the numerical results shown in Fig. [I] Note that, in 
one dimension, Eq. ^ describes a harmonic chain with 
random static friction [29[. The process described above 



is simply the transmission of forces when such a chain is 
pushed at its ends. 

As we now describe, the generalization of this argu- 
ment to two dimensions will account for the critical cur- 
rent of the bilayer, as measured by Tiemann et al. 0]. 
(We see similar behavior in two-dimensional simulations, 
albeit with large disorder fluctuations.) In our scenario 
of saturated domains, the static state only breaks down 
when the final domain in the sample exceeds threshold. 
Therefore, we may determine the critical current by set- 
ting / = 1 everywhere in Eq. ([5]) . Integrating over space, 
we see that the critical current, defined as the total in- 
jected current at threshold, functionly counts the number 
of domains in the sample, and is independent of the pre- 
cise geometry. We find, for its order of magnitude, 



ep s J 



(7) 



where S is the sample area in 2D and sample length in 
ID. Note that this is a natural form for an extensive crit- 
ical current, composed of the system size, the character- 
istic length Ld, and the microscopic current scale ep s /h. 
The same form can be seen in Eq. (13) of Ref. [2^. How- 
ever, it is not clear whether the result there applies to a 
bilayer driven at its boundary. 

For Eq. ([7]) to apply, the sample should comprise many 
domains. If a dimension L x is smaller than the domain 
size Ld then the second term in Eq. ^ should be multi- 
plied by \J Ld/ L x , because the total tunneling current of 
the domain is cut off at the sample width L x . This gives 

(quasi-lD: L, « L d « L y ) . (8) 

Ld L>d 

Similarly, for a sample containing only a single domain, 





(for L x ,L y < L d ). 



(9) 



To compare Eq. (J7J) with the experiments, we start 
from the microscopic theory [30l | for the zero-temperature 
values of the stiffness, p°, and order parameter, m x = 
(cos cf>), for a homogeneous bilayer. The latter renormal- 
izes the tunneling strength t° so that t° = Aom x /2irl 2 
where Ao is the single-particle tunnel splitting. This 
theory does not predict accurately the critical layer sep- 
aration for the loss of interlayer coherence but should 
be reasonable away from the critical point. Owing to 
the strong charge disorder, we expect the incompressible 
quantum Hall phase to occupy a small fraction of the 
sample, with the remainder occupied by puddles of com- 
pressible electron liquid 
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We suppose that the 
incompressible phase forms a network [18[ of channels of 
size I separating puddles of size £ ks dd, the distance to 
the dopants. Thus, the effective parameters in the energy 
[Eq. (TT])] for the disordered bilayer should be t ~ (l/dd)t° 
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and p s ~ (l/dd)Ps (~ 20mK at a layer separation d w i). 
The microscopic current scale is then ep s /h s=s 0.5 nA. 

Figure [5] shows our estimates for the domain size 
and the critical current I c , as functions of the ratio of 
inter layer separation to magnetic length, d/l. We see 
that the domains (L\ < 0.01 mm 2 ) are indeed not larger 
than the samples areas of (0.01 — 1) mm 2 so that the 
results are consistent with area scaling for I c . Moreover, 
I c has the correct order of magnitude compared to the 
observed values [§] of 0.1 — 10 nA. 

Interestingly, I c has a peak as a function of d/l which is 
also suggested in the experimental data [1] . This feature 
appears robust: it arises from the increase in Ld as d/l 
is reduced, caused by the increase in p s . However, the 
peak position depends on the precise variation of model 
parameters with d/l. For example, the variation in p s 
may cause some variation in vortex density [l7j and hence 
£, pushing the peak to smaller values of d/l. 




0.9 1.1 

d/l 



FIG. 2: Predicted critical current (solid, left axis) and domain 
size (dotted, right axis), for area S = 0.1mm 2 , d = 28 nm, 
£ (~ dd) — 200 nm, and bare tunneling Ao = 150 fiK. The 
stiffness and tunneling renormalization are taken from theory 
[30l |. and scaled by the network geometrical factor (see text). 

Our theory also has implications for the effect of an 
in-plane magnetic held Bu . This introduces a lengthscale 
in = hc/eB\\d, which is the length of a loop enclosing a 
flux quantum in the cross-section of a bilaycr. Within our 
theory, the system should be insensitive to Bu unless lii < 
£. At such fields, circulating tunneling currents are set up 
between the two layers, so that the net tunneling current 
in a domain is reduced. Therefore, coherent tunneling 
should be suppressed for B\\ > hc/ed£ ~ 0.7 T, which 
is consistent with experiments [3l[ where the enhanced 
tunneling decreases above 0.5 T. 

In conclusion, we have presented a theory of the crit- 
ical intcrlayer currents in a disordered quantum Hall bi- 
layer with static pinned vortices. We find that, because 
the current is injected at the boundary, coherent tunnel- 
ing saturates in the current-carrying region, leading to a 
Bean critical state. This results in an extensive critical 
current for sufficiently large samples (in contrast to the 
clean limit [HI where area scaling holds for small sam- 
ples) . The magnitude of the critical current is consistent 
with experiments. We predict that area scaling does not 
hold when the samples become smaller than the phase- 
pinned domains, and also that the critical current peaks 
in the interlayer coherent phase. 
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